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THE  ANGULAR  RESOLUTION  OF  MULTIPLE  TARGETS 


Ly 


J.  R.  Sklar 


F .  C .  Schweppe 


ABSTRACT 


The  angiilar  accuracy  chtalnahle  from  an  aperture  of  fixed  size 
Is  considered,  with  emphasis  on  the  multiple  target  case.  By  means 
of  the  Cramer -Rao  Inequality,  a  lower  hound  to  the  measurement 
variance  Is  computed.  In  addition,  the  degradation  due  to  uncer¬ 
tainties  In  the  number  of  targets  present  Is  considered.  Loss  of 
accuracy  Is  small  until  the  targets  approach  to  within  one  beamwldth, 
at  which  point  the  degradation  becomes  severe. 


This  technical  docvnnentary  report  is  approved  for  distribution. 


Air  Force  Lincoln  Laboratory  Office 


1-  Introduction 


Resolution  of  multiple  targets  has  always  been  of  some  concern  to  the  radar 
designer,  but  the  recent  interest  in  discrimination  has  given  a  new  importance 
to  be  problem.  In  a  multiple  target  environment,  energy  from  several  reflectors 
impinges  on  the  receiving  antenna  and  as  a  result,  the  data  processing  involves 
both  the  extraction  of  salient  features  of  the  received  signal  and  the  assign¬ 
ment  of  these  features  to  their  respective  targets .  Range  and  range  rate 
resolution  can  be  effected  by  the  design  of  the  radar  signeil  in  the  time  dimension 
and  angular  resolution  can  be  achieved  through  design  in  the  spacial  dimensions 
(antenna  design) .  However,  there  is  at  least  an  order  of  magnitude  difference 
between  the  angular  and  range  resolution  capabilities  of  today *s  more  sophis¬ 
ticated  radars.  The  possibility  that  this  discrepancy  is  at  least  partly  due 
to  a  failure  to  utilize  all  available  angular  information  in  the  return  signal 
is  indicated  by  the  success  of  the  monopulse  lobe  comparison  technique  for  beam 
splitting  when  only  a  single  target  is  present,  and  its  abrupt  failure  when 
there  are  more  than  one  in  the  same  range  resolution  interval.  Thus,  the  ques¬ 
tion  of  the  existence  of  a  middle  ground  in  which  a  moderate  degree  of  beam 
splitting  in  presence  of  multiple  targets  is  possible  arises . 

This  report  investigates  the  amount  of  angular  information  inherent  in  a 
received  signal  by  the  use  of  the  Cram^r-Rao  or  Information  Inequality.  Boxmds 
on  the  minimum  obtainable  angular  accuracies  are  derived  for  both  the  one-  and 
two-target  case.  The  detection  problem  of  deciding  the  number  of  targets  present 
is  also  discussed  and  limited  results  given. 

2*  Preliminaries 

The  following  two-dimensional  radar  model  is  used.  The  three-dimensional 
case  is  a  direct  generalization.  Several  targets,  say  Q,  are  reflecting  energy 
continuously  toward  a  linear  aperture  of  fixed  length,  X.  The  incident  plane 
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wave  is  sangjled  at  an  arbitrarily  large  number  of  points  on  the  aperture  and 
noise  is  added  to  the  samples.  Since  the  transmission  mediimi  is  linear,  the 
contributions  to  the  sample  value  from  each  target  can  be  summed,  and  therefore 
the  samples  are  a  2  component  vector  consisting  of  the  quadrature  components 

Q 

z^(x)  Aj  sin  (cpj  +  <I>(x,aj))  +  w^(x)  0  <  x  SX 

J=1 


Q 


>(x)  cos  (cpj  +  <I>(x,aj))  +  WgCx) 

J=1 


(1) 

th 


where  A.  and  cp  .  are  the  magnitude  and  phase  of  the  radiation  from  the  j 
J  J 

target,  CX  is  its  angular  position,  x  is  the  location  of  the  sample  point 
J 

within  the  aperture,  and  <I>  is  a  fimction  of  x  and  0!  which  measiires  the  phase 

J 

difference  between  sample  points  due  to  the  geometry.  (See  Fig.  1.)  If  X 
is  the  wavelength  of  incident  plane  wave,  then 

cos 


Wj^(x)  and  w^Cx)  are  Gaussiein  random  variables  of  zero  mean  representing  the 
noise.  Since  w^(x)  and  w^Cx)  are  orthogonal  components  of  the  rf  noise,  they 
are  independent.  In  addition 


E  [w^(x)w^(y)]  = 


X  x'y 

X  =  y 


i  =  1,2 


The  anqjlitudes.  A.,  the  phases  cp.,  and  the  angles  cc ,  are  considered  to  be 
J  J  0 

unknown  parameters;  i.e.,  no  a  priori  distribution  for  them  is  assumed  known. 
Thus,  there  are  3Q  unknown  parameters  and  the  problem  is  of  the  basic  form 
Zl(x)  =  ,e^^)  +  w^(x) 


^2^^^  "  +  WgCx) 


(2) 
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TARGET  j 


3-28-4804 


X 

-• -  LINEAR  APERTURE  OF  LENGTH  X - 

Fig.  1.  Target  Geometry. 


where  the  € ' s  represent  the  unknown  parameters ,  and  the  f .  *  s  are  known 

J 

functions  of  x  and  the  e’s.  We  will  define  ^j(j  ==  as  the  estimate 

of  the  obtained  from  the  observations  z^(x)  ,  Z2(x)  ;  0  <  x  ^  X. 

We  consider  the  case  where  the  are  observed  at  N  points  of  the 

aperture,  equally  spaced  at  intervals  6  =  X/N. 

Consider 


p(z(i) , . . . ,z(n) 1 

th 

where  z(k)  is  a  two-dimensional  vector  representing  the  k  sample,  and  p  is 

the  probability  of  the  observed  sample  conditional  on  the  true  values  of  the 

unknown  parameters.  If  one  desires  to  obtain  \inbiased  estimates  of  a  set  of 

unknown  quantities  based  on  observing  a  set  of  random  variables  z(k) ,  one 

J 

can  place  a  lower  bound  on  the  variances  of  the  estimators  €  by  means  of  the 

J 

Craxner-Rao  or  Information  Inequality.  Two  of  many  derivations  of  this 
fundamental  result  are  foiand  in  Ref.  2,  Chap.  32  and  Ref.  3j  Chap.  12  (the 
terms  Carmer-Rao  or  Information  Inequality  are  not  employed;  see  instead  the 
discussions  on  efficiency).  In  the  case  of  multiple  parameter  estimation, 
this  inequality  can  best  be  stated  in  matrix  notation.  Let 


I.  .  =  -E 

ij 


^  log  p  (_z(  l) , . .  .^(N)  )  j  e^,...e2Q 

^  J 


be  elements  of  a  matrix  [l] ,  and  let 


"kj  =  -  'i*  ('j  -  "j)] 

be  the  elements  of  a  matrix  .  Then  the  Cramer-Rao  Inequality  states  that 

/]  - 

is  positive  definite,  which  implies  that 


-k- 


(3) 


This  result  is  derived  from  the  requirement  that  the  estimates  be  lanbiased  and 
a  generalized  Schwartz  inequality  and  is  subject  to  the  condition  (usually 
present  in  practice)  that  [l]  be  positive  definite.  Thus  the  desired  lower 
boiind  can  be  obtained  by  calculating  the  necessary  derivatives  and  inverting 
the  resulting  matrix. 

In  the  sequel  it  will  be  useful  to  employ  the  following  alternate  inter¬ 
pretation  of  [l]  Consider  Eq.  (2).  Expand  the  f in  a  Taylor 

series  about  some  set  of  e’s,  €?;...  Then 

1 


Az^(x)  = 


A€j  +  W^(x)  +  ji  =  1;2  (4) 


J=1 


where  the  partial  derivative  is  evaluated  at  the  €°'s.  is  the  remainder. 


and 


Az^(x)  =  z^(x)  -  f^(x,€°,...,€°  ) 


If  the  are  close  enough  to  the  true  values  to  allow  the  remainder, 

to  be  neglected.  Eq.  (4)  can  be  considered  as  a  linear  regression  problem  with 

unknown  coefficients,  Ac..  A  minimim  variance,  lanbiased  estimate  (or  a  maximum 

J 

likelihood  estimate)  for  this  problem  has  a  covariance  matrix  equal  to  [I] 
[References  2  and  3  give  the  linear  regression  equations  which,  when  applied 
to  Eq.  (4),  give  [l]  ^  for  the  Gaussian  case.]  Therefore,  [l]  ^  can  be  con¬ 
sidered  as  either  a  bounding  covariance  matrix  or  the  result  of  a  linearized 
or  small  error,  error  analysis.  Thus,  for  signal-to-noise  ratios  which  are 
large  enough  to  Imply  small  errors,  the  inequality  of  Eq.  (3)  becomes  equality. 
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single  Target 


The  first  case  to  be  considered  is  a  single  target  at  unknown  angle  a  with 
unknown  anplitude  A  and  phase  0.  The  target  is  known  to  exist.  When  the  noise 
is  gaussian,  with  variance  and  zero  mean,  the  observation  of  N  samples  equally 

X 

spaced  at  6  =  —  along  the  aperture  leads  to  a  conditional  probability  for  the 
sacQjles  of 


p(z(6),z(26),...  |a,A,0) 


N  r- 


-z 

i-1 

N 


z  ( ib)  -A  sin(0  +  ~  cos  a) 

X  \ 


2c 


“  ^j^2(i6)  -A  cos(0  + 


2TTi6 


cos 


i=l 


from  Eq.  (l) . 

After  taking  the  proper  derivatives  and  averages  we  obtain,  for  N  large, 

,2 


rii  = 


N  A^  Z' 2nx'f 

5 


.  2 

sin  a 


N  A  /  2ttx\  ,  ^ 

-2  tI,— 


W 

CT 


N  A^  /  2txX'N  .  „ 

2  tV— 


NA 


in  which  the  indices  1,  2  and  5  correspond  to  a,  A  and  0  respectively.  From 
this  we  calculate 
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12  6 


Consequently,  the  variance  in  the  angle  measurement  is  lower  bounded  by 


2 


a 


2 

sin  a 


(5) 


The  fractional  bearawidth  error  is  therefore  bounded  by 


m 


where  R  =  is  the  signal -to -noise  ratio  of  the  same  aperture 

element  outputs  are  added  and  BW  =  is  the  bearawidth  for  the 
'2^^ 


measured  to  the  (  —  relative  power  level.  We  see  that  this  is 

\TT 

that  obtained  by  monopulse  lobe  comparison  techniques  in  Ref.  1. 


if  all  array 
aperture 
comparable  to 


h .  Two  Targets 

The  same  procedure  can  be  followed  in  the  two -target  case  when  exactly 
two  targets  are  known  to  exist.  The  conditional  probability  expression  is 
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P(z(6),z(26),... 


(2tt) 


N  2N  )  “  ? 

^  /  i=l 


z^(i6)  ~  sin(<i^  +  -^16)  cos  a^) 


-Ag  sinCOg  +  ^i6)  cos  a^) 


N 

i 


^2^  iB) 


i«l^ 


-  cos(o^  +  ^iB)  cos  a^)  -  Ag  8ini(4>2  +  *^2^ 


Taking  derivatives  and  computing  averages  the  result  is 

I  -  [D]  rrj  p]  (6) 

where  [P]  is  a  6  x  6  matrix  given  in  Fig.  2  and  p]  is  the  6x6  diagonal 
matrix, 


2t«  .  ,  ^ 

-  A.  sin  a. 

X  1  1 


m- V 


2ttX  .  .  ^ 

—  A^  sin 


-6 


0 


Cn  = 


V5 

^21 

0 

1/2 

^61 


21 

1/5 

■^41 

0 

^6l 

1/2 


■^41 

1 

^4l 


^65 


"41 


45 


"65 


l/2 


'6l 


'65 


'45 


^6l 

1/2 

-a 


65 


'45 


21 


(aX) 


^  cos(aX  +  b)  +  7-4y  sln(aX  +  b) - sln(aX  +  b)  +  — ^  sin  b 

(aX)5 


(aX)- 


^),i  “  r^T  cos(aX  +  b) - sin(aX  +  b)  +  — sin  b 

(axr  (aX)^ 

“  ZiVN  sin(aX  +  b)  +  — cos(aX  +  b) - cos  b 

(axr  (axr 

“1,3  ■  =1”  ““ 

‘63  ‘  +  >>)  -  ““ 

2tt  ,  . 

I  =  —  (cos  —  cos  ct^) 


=  «?!  -  «P2 


Figure  2  -28-4805 

[n  Indices  1,  2,  5,  4,  5,  and  6 
correspond  to  parameters  a^,  A^, 

Ag,  (p^,  and 
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This  factorization  of  the  fll  matrix  is  attractive  since  it  places  in 
evidence  the  dependence  of  the  elements  of  [l]”^  on  the  various  parameters. 
(Recall  the  inverse  of  a  diagonal  matrix  is  a  diagonal  matrix  with  the  elements 
inverted  and  also  [ll  ^  =  [D]  ^  [r]  ^  [D]  ^.)  Thus  we  may  note  that  the 
element  of  [I^  ^  corresponding  to  the  variance  of  aii5)lltude  is  independent  of 
the  an5)lltude.  This  fact  will  be  important  in  the  detection  discussion  of 
Section 

Computing  the  inverse  of  the  fll  matrix  is  difficult,  but  when  we  are  only 
concerned  with  the  angular  accuracy,  only  part  of  the  inverted  matrix  is  re¬ 
quired.  This  can  be  obtained  with  a  somewhat  less  tedious  calculation  by 
partitioning  the  [r]  matrix.  Althou^  it  is  then  possible  to  obtain  the 
variance  of  the  angular  measurement  in  closed  form,  the  expression  is  too 
lengthy  to  be  very  useful.  Therefore,  the  angular  variance  is  plotted  as  a 
function  of  the  target  separation  in  Fig.  5  with  the  phase  difference  b  =  <t>^ 

as  a  parameter.  From  this  normalized  curve  results  for  various  values  of  a 
signal -to -noise  ratio  and  angular  position  can  be  obtained. 

In  addition,  the  ratio  of  the  variance  in  the  two -target  case  to  the  var¬ 
iance  in  the  one -target  case  is  also  plotted.  One  can  observe  that  the  results 
are  essentially  identical  for  target  separations  greater  them  one  beamwidth 
(BW  -  f  that  for  smaller  separations  the  measurement  error  in  the  two- 
target  case  becomes  large  very  rapidly  as  the  separation  decreases.  The  phase 
difference  b  causes  a  wide  variation  in  the  form  of  the  curves.  However,  for 
all  values  of  this  parameter,  the  variance  is  radically  larger  if  the  separation 
is  less  than  about  one  beamwidth.  Consequently,  any  attempt  at  beam  splitting 
for  separations  much  smaller  than  one  beamwidth  will  prove  futile. 

3.  Detection 

The  preceding  two  sections  have  considered  estimating  accviracy  for  multiple 
targets  under  the  assuinptlon  that  the  exact  number  of  targets  present  is  known. 
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Fig.  3.  Angular  Accuracy  in  Presence  of  Two  Targets. 


Detection  is  the  problem  of  deciding  the  nmber  of  objects  present.  Certain 
general  aspects  of  multiple  target  detection  are  briefly  reviewed,  and  a 
measure  provided  on  the  relative  detection  capability  in  the  one-  and  two- 


target  case. 

A  standard  approach  to  the  detection  problem  employs  the  theory  of 
hypothesis  testing.*  In  the  case  where  at  most  one  target  may  exist,  it  is 


hypothesis) .  The  alternate 


hypothesized  that 


hypothesis  is  that  a  target  does  exist  (H^  is  the  alternate  hypothesis) .  If 
A  is  the  amplitude  of  the  signal,  these  hypotheses  are  equivalent  to 


A  =  0 
H,:  A  |^  0 


If  up  to  tvo  targets  are  possible >  it  is  then  necessary  to  check  the  possi¬ 
bilities  of  zero,  one  and  two  targets.  Such  a  multi-level  decision  can  be 
built  out  of  two-level  tests  such  as  used  for  the  one-target  problem.  For 
example,  a  first  test  might  hypothesize  that  no  targets  are  present  vs  the 
alternate  hypothesis  that  at  least  one  target  is  present.  That  is 


At  least  one  ^  0  j  =  1,2 


(7) 


targets  are  present  is  rejected,  one  could  then  hypothesize  the  existence  of 
exactly  one  target  vs  the  alternate  that  exactly  two  targets  exist.  That  is 


Hq-.  0  ,  Ag  =  0 

A^  0  ,  0 


(8) 


*The  following  discussion  is  a  very  special  application  of  the  general 
statistical  theory.  References  2  and  3  ai'e  two  of  many  which  contain  far 
more  extensive  discussions. 
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There  are,  of  course,  many  other  possible  sequences,  and  the  above  is  only  an 
example . 

The  basic  problem  is  thus  the  testing  of  some  hypothesis  Hq  vs  some 

alternate  The  test  is  based  on  some  test  statistics,  where  ^  is 

some  function  of  the  observed  data.  There  are  many  possible  test  statistics, 

two  of  which  have  particularly  nice  properties.  The  first  of  these  is  based 

directly  on  the  estimated  amplitudes  (X,)  of  the  signals.  For  example,  in 

J 

the  one-target  case,  the  hypothesis  that  no  target  exists  (Hq)  is  rejected 
if  the  magnitude  of  the  estimated  amplitude  exceeds  some  chosen  value.  In 
the  multiple-target  case,  similar  tests  are  enrployed  on  the  vector  of  the 
estimated  amplitudes.  The  second  test  statistic  is  the  likelihood  ratio; 
i.e.,  the  ratio  of  the  maximum  likelihood  attainable  under  to  the  maximum 
likelihood  attainable  -under  H^.  (See,  for  example.  Ref.  3»)  The  amplitude 
statistic  is  more  powerful,  but  the  likelihood  ratio  test  is  independent  of 
the  signal- to-noise  ratio.  We  shall  confine  discussions  -to  the  amplitude 
statistic. 

To  be  explicit,  assume  the  signal- to -noise  ratio  is  high  enough  so  that 
the  variance  given  by  the  Cramer-Rao  Ineqxiality  can  actually  be  realized  or, 
equivalently,  that  we  are  dealing  -with  thesystem  of  linear  eqxiations  given  by 

/ 1  \  r  ^  1 

Eq.  (4).*  The  vector  estimate  IA^,A2J  is  then  a  -two-dimensional  Gaussian 
random  variable  wi-th  covariance  matrix  given  by  -the  corresponding  elements 
of  the  matrices  of  Let  denote  this  2x2  sub-matrix  of 

If  it  is  desired  to  test  the  hypothesis  of  Eq.  (7)^  the  test  statistic  ^  is, 
in  matrix  notation,  i — 


t"  = 


/V 

A 

■ 

I 

1 

2. 

-1 

AA 


A 

A, 


A 

A. 


*This  approach  is  similar  to  the  investigation  of  the  asyniptotic  beha-vlor  of 
test  statistics  els  enployed  in  the  statistical  literature. 
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2 

The  equation,  ^  =  constant,  defines  an  e3J.ipse  of  constant  probability  in 

A  A  2  2 

Ag  space.  Thus  the  hypothesis  test,  which  rejects  when  5  ^ 

the  hypothesis  is  rejected  whenlA^,  A^llies  outside  the  constant  probability 

2  A 

ellipse  specified  by  .  Let  denote  the  variance  of  A^  as  evaluated, 
using  the  corresponding  diagonal  term  of  For  the  case  of  Eq.  (8),  the 

test  statistic  is 


which  is  the  one -dimensional  version  of  the  preceding  case. 

The  question  of  prime  interest  in  this  paper  is  how  well  detection  can  be 
performed.  Unfortunately,  a  complete  emalysis  reqtiires  the  choice  of  explicit 
sequences  of  tests,  and  this  choice  depends  on  the  particular  problem  of 
interest.  In  addition,  the  calculation  of  the  resulting  probabilities  of 
the  vario\is  types  of  error  is  somewhat  laborloios.  Thus  we  shall  merely 
indicate  the  relative  degradation  in  detection  capability  which  results  from 
the  presence  of  a  second  target.  Consider  the  case  of  Eq.  (8) ,  •vdiere  the 
existence  of  one  target  is  definitely  known  (A^  ^  O) ,  and  the  possible  exis- 
tence  of  a  second  is  to  be  tested.  The  test  statistic,  g  ,  is  ^ 


5  •» 


where  q  is  a  Gaussian,  zero-meem  random  variable  of  laiity  variance.  Now,  con¬ 
sider  the  problem  of  deciding  whether  there  are  one  or  no  targets.  Ebcactly  the 
same  test  is  used,  except  that  is  replaced  by  the  obtained  from  the 
corresponding  element  of  [l]  as  evalmted  in  Sec.  3^  Fq.  (5).  Thus,  for  a 
fixed  false  detection  probability,  the  probability  of  detecting  a  second  target 
of  anplitiide  Ag  in  the  presence  of  another  known  to  be  present,  eqmls  the 
probability  of  detecting  a  single  target  of  amplitude  A  when 


The  ratio,  therefore  measures  the  degradation  of  detection  performance 

caused  by  the  presence  of  a  second  target. 

As  with  the  case  of  the  angular  accxiracies  of  Sec.  4,  the  closed  form 

2  2 

solutions  are  too  complex  to  be  very  useful.  Thiis  /a^  is  plotted  in  Fig.  4 
as  a  function  of  target  separation  with  the  phase  difference  b  as  a  parameter. 

6.  Discxission 

Figures  3  and  4  constitute  the  major  resvilts  of  this  study  and  show  that, 
as  target  separations  increase  beyond  one  beamwidth,  the  information  loss 
caused  by  multiple  targets  rapidly  disappeeurs.  Since  these  curves  resulted 
from  a  limited  analysis  on  a  siniplified  model,  we  must  discuss  their  limita¬ 
tions  before  drawing  actual  conclusions. 

We  considered  only  errors  due  to  additive  noise.  Thus  the  antenna  pattern 
is  assumed  to  be  known  exactly.  This  assumption  results  in  the  angular  accvuracy 
of  one  target  being  independent  of  the  cross  sections  (anip3JLtude)  of  other 
targets.  With  an  exactly  known  antenna  pattern,  this  is  not  unreasonable^  since 
the  effect  of  other  targets  can  be  "subtracted  out."  In  practice,  however, 
uncertainties  in  antenna  pattern  can  caixse  a  large  cross  section  target  to 
sizeably  degrade  estimates  from  a  nearby  smaller  target.  In  addition,  an 
ineqviality  was  ixsed  in  place  of  an  exact  error  analysis.  We  can  hope  to 
approach  the  resulting  bound  only  with  a  sufficiently  high  signal -to-noise 
ratio.  Thus  Figs.  3  and  4  indicate  actual  capabilities  only  for  S/N  ratios 
high  enough  to  approach  the  bound  but  low  enough  to  keep  antenna  uncertainties 
from  predominating.  The  range  of  such  s/N  ratios  depends  on  the  particular 
problem  xmder  investigation.  However,  it  is  not  unreasonable  to  expect  some 
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Fig.  4.  Detection  Degradation  in  Presence  of  Two  Targets. 


systems  for  which  it  is  nonempty.  For  example,  the  same  basic  assumptions 
hold  in  both  the  one-target  and  two-target  cases  and  the  single  target  bound 
is  often  approached  in  practice. 

We  applied  the  basic  theory  to  a  very  simple  model:  a  one-dimensional 

aperture  with  two  stationary  targets  and  white  noise.  The  two-dimensional 

aperture  and  more  than  two  targets  merely  introduces  more  parameters  into 

the  problem.  This  is  also  true  for  moving  targets  as  the  unknown  phase  0 

til 

of  the  return  signal  from  the  j  target  becomes  a  parameterized  f\mction  of 
time  such  as 

iiij(t)  =  +  jiijgt 

to  acco\mt  for  range  and  range  rate  and  the  angular  position,  a  ,  becomes  a 

J 

similar  parameterized  function  of  time  to  acco\mt  for  the  angular  motion. 

A  moving  target  also  requires  specification  of  the  carrier  modulation.  The 
introduction  of  these  parameters  into  the  model  is  straightforward  and  would 
not  change  the  basic  method  of  analysis.  However,  it  would  greatly  enhance 
the  conputational  difficulties  associated  with  explicit  results.  The  simple 
model  is  valuable  as  it  indicates  the  target  separations  for  which  the  angular 
resolution  problem  can  be  ignored  and  either  single  target  analysis  used  or 
range  and  range  rate  mioltiple  target  resolution  studied  by  itself.  Such  a 
dichotomy  provides  a  sinplified  but  very  useful  picture  of  the  over-all  multiple 
target  problem.  The  incorporation  of  n  order  Markovian  observation  noise 
changes  the  eq\iation  for  information  to  a  nonlinear  differential  eqmtion,  which 
can  be  nvnnerically  integrated  on  a  digital  conputer.  However,  our  basic  problem 
is  not  well  enough  defined  to  justify  a  nonwhite  noise  model. 

The  most  inportant  limitation  in  the  scope  of  the  study  is  the  fact  that 
we  have  not  considered  an  actual  data  processing  system.  If  conditions  are 
such  that  the  bound  can  be  approached,  the  classical  maxlmimi  likelihood 
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estimation  technique  provides  the  desired  estimates.  If  we  are  not  worried 
about  computation  time  and  cost,  such  estimates  can  be  obtained  vising  a  large 
digital  computer.  However,  for  most  situations,  such  a  procedure  is  grossly 
impractical,  as  real  time  estimates  are  required  and  these  imply  primarily 
analog  data  processing.  Unfortunately,  a  straightforward  einalog  implementation 
of  the  maximum  likelihood  technique  for  several  targets  requires  a  large 
amoung  of  hardware,  and  the  various  avenues  which  might  lead  to  a  practical 
analog  implementation  have  not  yet  been  explored.  Information  is  of  no  value 
if  it  cannot  be  extracted. 

Let  us  now  summarize  what  Figs.  3  and  4  actually  ing)!^.  They  provide  a 
bound  on  possible  system  performance.  More  important,  there  is  good  reason  to 
believe  that  some  systems  do  or  will  exist  for  •vdiich  this  performance  can  be 
approached.  Thus,  for  more  than  a  beamwidth  separation,  multiple-target 
capability  close  to  that  of  single  target  monopulse  lobe  compression  is  a 
distinct  possibility,  provided  a  practical  data  processing  procedure  can  be 
developed  and  implemented.  A  prime  puipnase  of  this  report  is  to  motivate 
investigations  on  such  data  processing  techniques. 
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